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' We consider random energy landscapes constructed from d-dimensional lattices or trees. The 

»' ' distribution of the number of local minima in such landscapes follows a large deviation principle 

and we derive the associated law exactly for dimension 1. Also of interest is the probability of the 
^ ' maximum possible number of minima; this probability scales exponentially with the number of sites. 

Q I We calculate analytically the corresponding exponent for the Cayley tree and the two-leg ladder; 

for 2 to 5 dimensional hypercubic lattices, we compute the exponent numerically and compare to 
the Cayley tree case. 
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' For many materials that are glassy [H local minima of the energy (or of the free energy) trap the system for long 
^ , times, leading to subtle equilibrium and out of equilibrium properties. Energy landscapes provide a simple conceptual 

a framework for modeling these systems but in fact their use goes much beyond that. For instance the vacua of string 
I ] theories are expected to proliferate enormously, and the problem of estimating the number of local minima [3, [s] 
1 of string energy landscapes is still an open problem. Other examples include rugged (random) energy landscapes 
I in evolutionary biology Q , quantum cosmology [5j , manifolds in random media Q , glassy systems ^|, 0] random 
^ , potentials [1, Q and the associated problems in random matrix theory [3, [13, [HI ■ Typically these systems consider a 

■ particle, a configuration of particles, or even a manifold, subject to a random potential. One extreme case for its visual 

a simplicity is that of a point particle in a random Gaussian potential; on the opposite extreme are energy landscapes 
I ' associated with the configuration of a many body system. Examples in this last category include: (1) the p-spin glass 
i-rt ! model |12|. which in the limit of large p reduces to the far simpler random energy model of Derrida [III; (2) atomic 

■ clusters|14| and other glassy systems with no quenched disorder fis'l ; (3) random manifolds in random media Q . 
O For any landscape, it is desirable to know the statistical properties of the minima (or more generally of the saddle 

, points). A quantity frequently considered is the expected number of minima as a function of their energy It r nay 

also be of interest to consider how the set of minima are organized topologically, e.g., whether the barrier tree [Iq 
. is ultrametric. Our focus here is to better understand the statistics of the total number of local minima in random 
^ ■ energy landscapes; our underlying space is a regular (Euclidean) lattice on each site of which resides a random energy. 
[ Let M denote the total number of local minima for given values of the energies on each site. Evidently M is a 

■ random variable as it varies from one realization of the landscape to another. We are interested in the statistics of 
Q,^ ' M as a function of the number N of lattice sites. Similar questions were studied recently in the context of random 

i permutations ^17i,il8j, ballistic deposition il9i |, and in simple models of glasses j20j . In this paper, we provide a number 
of analytical results on the distribution of the total number of local minima for random energy landscapes on several 
lattices; the moments of M are easily derived, so our focus concerns mainly the probability of large deviations of M 
from its mean value, i.e., the probabilities of atypical configurations. Furthermore we ask what is the probability of 
M being at its maximum, corresponding to the limit where the minima are maximally packed on the lattice. 
^ • The paper is organized as follows. In Sect. |T] we specify the model and show that the statistical properties of the 
minima are independent of the individual distribution of energies as long as these energies are independent from site 
to site and are drawn from a continuous distribution. In Sect. |TT] we cover some of the simplest properties of the 
statistics of M and formulate the large deviation principle. Then we derive the closed form expression for the large 
Q , deviation function in the case of the one dimensional lattice in Sect. IIIII The focus of the rest of the paper is the 
maximum packing problem. In Sect. II Vl we determine the probability of maximum packings for several solvable cases, 
namely Cayley trees and a two-leg ladder. The case of d-dimensional hypercubic lattices is then treated by numerical 
computation in Sect.|Vl Finally, these different results are discussed in Sect. I VII and some closing remarks are given. 



I. THE MODEL 



We start with a regular lattice; on each site i lives a random energy Ei, drawn independently from site to site from 
a common distribution p{E). We assume p{E) is continuous and normalized to unity, i.e. p{E)dE = 1. For any 
choice of the set of E'i's we obtain an energy landscape, i.e., a topological space (defined through nearest neighbors 
on the lattice) with an energy for each element. For a given realization of the landscape, a site i is a local minimum 
if Ei < Ej where j denotes any nearest neighbour sites of i. Hereafter we shall denote by M the number of local 
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minima on this landscape. 

For any given realization of the landscape, one can formally express M as 



N 



^'^^Y. n (^iE,-E,) (1) 

^=1 3l<i]> 

where the product runs over all nearest neighbours j of site i and Q{x) is the Heaviside theta function. We are 
interested in computing the probability distribution P{M, N) of the total number of local minima (or equivalently 
that of the maxima) . 

The first important observation is that the distribution P{M, N) is completely independent of the energy distribution 
p{E). To see this, we formally express the distribution as 



P{M,N) = f ... f sIm-J^HHEj-eA Y[piEk)dEk 

\ 1=1 j J k=i 



(2) 



where 5{x) is the Dhac delta function. Next we make a change of variable for each k 



Xk = / p{E)dE (3) 

J — oo 

Clearly Xk is a monotonically increasing function of Ek. Therefore, 9 {Ej — Ei) = d{xj — Xi). Moreover, since p{E) is 
normalized to unity, the variable x varies from to 1 and so Eq. ^ simply becomes 

(4) 



P{M,N)^[ dxi... [ dxNS\M^y2T\0{xj-x,) 
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We see that the energy distribution p{E) simply drops out and P{M,N), for arbitrary p{E), is universal and is the 
same as when the 'new' energy variable Xi is drawn independently from a uniform distribution over x G [0, 1]. Thus, 
the model is simplified: at each lattice site lives a random number Xi € [0, 1] drawn independently from a uniform 
distribution. We will refer to this model as the random minima model, for which we want to compute the distribution 
P{M, N) of the total number of local minima. 

Following the mapping to the uniform distribution, it follows that the random minima model is just the continuous 
version of the permutation generated landscape recently studied by Hivert et. al. (loj . They considered a set of 
integers {1, 2,3,..., N}. Each of the A^! permutations of this set defines a random energy landscape and they occur 
with equal probability. Thus the energy at a site is now a discrete integer drawn uniformly from the set [1, 2, 3, . . . , iV] 
with equal probability Thus the statistics of the minima in the permutation generated landscape will be identical 
to that of the random minima model. 



II. GENERAL PROPERTIES OF THE DISTRIBUTION OF M 



The mean and the variance of M are relatively straightforward to compute on an arbitrary lattice since they depend 
only on the local properties of the landscape. To see this, let us define the variable 



(5) 



where j runs over the nearest neighbours of i and the x^'s are independent random numbers in [0, 1] drawn from the 
uniform distribution. This rji is an indicator function which is 1 if site i is a local minimum and otherwise. Then it 
follows that 



N 



(6) 



Taking the average in Eq. ^ and using the translational invariance of the lattice, it follows that 
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where n is the number of nearest neighbours of any site, i.e., the co-ordination number of the lattice. This result 
also follows trivially from a combinatorial argument on the permutation landscape: consider the site i with its n 
neighbours. The number of ways one can arrange the integers on these (n + 1) sites with the restriction that z is a 
local minimum is nl clearly. On the other hand, the total number of unrestricted configurations is [n + 1)! so the 
probability that the site i is a minimum is simply nl/(n + 1)! = l/{n + 1). 

The calculation of the variance takes a few more steps. Squaring Eq. ([6]) and taking the average gives 



(M2)=^(,y,r;,). (8) 
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Note that the rji^s are correlated random variables but only over a short range. The calculation of the correlation 
function {rjirjj) can therefore be performed by hand (it only involves the calculation of simple integrals involving a 
maximum of 2n sites) . For example, in 1-d, one gets 

(7^ ^ (A/2) - {Mf as N ^oo. (9) 

This result has been obtained in various contexts before, such as in the calculation of the number of metastable states 
in 1-d Ising spin glass [IH and also in the context of the occupation time of a non-Markovian sequence in 1-d [2^ ]. 
Recently, the variance for the 2 — d square lattice was computed in the permutation model by Hivert et. al. [l^. For 
large TV, one gets 

In general, on any arbitrary lattice, for large N 

(M) = aN 

(M^) - {Mf = bN (11) 

where a and h are lattice dependent numbers that can be computed either using integrals or by combinatorics in the 
permutation model. 

Near the mean (M) and within a region \M — (M)| = 0{^/N), the distribution P{M,N) is expected to be a 
Gaussian |T§| with mean and variance given in Eq. (jlip 

P{M, N) ^ -=L= exp [-(M ~ aNf/2bN] . (12) 
vzTToJy 

However, for M far from the mean, one would expect deviations in P{M, N) from the Gaussian form. In this paper, 
we are interested in the probabilities of such large deviations, i.e., we are interested in computing the probabilities of 
occurrences of configurations that are far from typical. On general grounds, in the limit M — > oo, N oo but with 
their ratio M/N fixed but arbitrary, one expects that P{M, N) has the form 



P{M, N) ~ exp 



N 



(13) 



where $(?/) is a large deviation function. Now, on a given lattice, Af can take values from to a maximal number 
A/max < N . The upper limit follows from the fact that once a site is a local minimum, none of its nearest neighbours 
can be a local minimum. Thus there is a nearest neighbour exclusion principle for the minima. This constraint 
indicates that one can not pack arbitrarily large number of minima on the lattice and there is an upper bound on 
the number of minima. For example, on a bipartite lattice, one can pack at most A/max — N/2 local minima, one at 
every alternate site. In particular, on a square lattice the minima will be placed on a checkerboard pattern. Thus, in 
this case, the large deviation function $(1/) is defined for < y < 2/max = 1/2. On the other hand, on a Cayley tree 
with fi number of branches and N sites, we will see later that Mj^ax — I^N/ (fj, + 1), thus ?/max — fJ-/ifJ- + 1)- 

As mentioned above, in the vicinity of its mean, i.e., for \M — {M)\ = 0{Vn), the distribution P{M, N) is Gaussian 
as in Eq. (|12p . This indicates that the large deviation function $(j/) is quadratic near y = a 

'^(y) - ^ (14) 

such that P{M,N) w exp [—N{y — a)"^ /2b] ~ exp [—{M — aN^ /2bN] has the required Gaussian form. However, 
far from the mean, P{M,N) will have non-Gaussian tails, indicating a departure of ^{y) from the simple quadratic 
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form. A knowledge of the function ^{y) would then allow one to compute the probabilities of occurrences of atypical 
configurations, such as the probability of a configuration with the lowest number of minima (e.g. AI — 1) or the ones 
with the maximal number of minima {M = Mmax)- 

A particular focus of this paper will be to compute the probability of a maximally packed configuration, i.e., when 
M = Mmax- It follows from the general form in Eq. (fT3|) that the probability of this maximal packing should decay 
exponentially with the system size for large N, 

P{M^... N) ~ exp [-iV$(2/max)] ~ 7""^ (15) 

where 7 = exp [$(j/max)] is a lattice dependent constant. Of course, this expression also gives the scaling of the 
fraction of permutations which have Mmax local minima in the permutation landscape. It is then natural to interpret 
In 7 as the entropy cost per site to achieve maximally packed minima. 

We shall compute the constant 7 exactly for a number of lattices. For the 1 — d chain, we will show that 

7^1 = 1.57079... (16) 

Note that this result appeared before in the context of metastable states in 1 — d Ising spin glasses at zero tempera- 
ture Furthermore, the same (7r/2)~^ decay also appeared as the persistence probability of a 1 — d non-Markovian 
sequence [l^, H^l- For the 1 — d case, one can also compute the full large deviation function <&(?;) (see Section Hill . 
Another lattice where we can calculate 7 exactly is the Cayley tree with /i branches where we show that 

where B{x,y) is the Beta function. As expected, for ^ = 1, it reduces to the 1-d result in Eq. (flB]) . As fi increases, 7 
increases slowly and as/^— >oo,7^2. A Cayley tree with an infinite number of branches corresponds to a hypercubic 
lattice in d-dimensions in the d 00 limit. Thus, based on these two limiting results, one expects the constant 7 on 
any d-dimensional lattice to satisfy the bounds 

7r/2<7<2. (18) 

Our numerical simulations for 2 < d < 5 are consistent with these bounds. We have also been able to compute 7 
exactly for a two-leg ladder where we show that 



9a' 



1/3 

= 1.57657... (19) 

where a = 1.86635 ... is the smallest positive root of the Bessel function J_i/3(z) = 0. The result in Eq. is, of 
course, consistent with the general bounds in Eq. 
We now consider these successive cases in detail. 



III. EXACT LARGE DEVIATION FUNCTION IN ONE DIMENSION 



In a 1-d chain of size N, the distribution P{AI, N) of the total number of local minima was computed exactly by 
Derrida and Gardner [2l[, although they did not calculate the large deviation function ^{y) explicitly. However, from 
their result for the generating function, it is easy to derive the large deviation function by a Legendre transform. For 
the sake of completeness, we provide here a brief derivation of the 1-d result, albeit by a slightly different method. 

It is useful to define the generating function or the partition function Z{z, N) = J2m N)z^'^ where z is the 

fugacity or the weight associated with each minimum. For simplicity we consider an open chain of size N and let 
xn = a; be the value of the random variable at the 7V-th site. To write a recursion relation for the partition function, 
it is convenient to define two restricted partition functions: Zi{x, z, N) and Zo{x, z, N) denoting respectively the 
partition functions conditioned on the fact that the 7V-th site has value xn — x and that it is respectively a local 
minimum (i.e., xn < xn-i) or a local maximum {xn > xn-i)- Note that since we arc considering an open chain, 
the last site (N) has only one neighbour to its left, namely the {N — l)-th site. Knowing the restricted partition 
functions, one can compute the full partition function from the relation 



Z{z,N)= [ dx[Zi{x,z,N) + Zoix,z,N)] 
Jo 



(20) 
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The restricted partition functions satisfy a pair of simple recursion relations 

Z^{x,z,N) = z f Zo{y,z,N-l)dy+ j Zi{y, z, N - 1) dy (21) 
Zo{x,z,N) = f [Zo{y,z,N^l) + Zi{y,z,N~l)]dy. (22) 



These recursion relations can be easily understood by considering all possibilities when one adds a new site to the 
chain. If the new site is a minimum, we need to attach a factor z. On the other hand, if the (TV — l)-th site was a 
minimum and it ceases to be a minimum after the addition of the A'^-th site, we have to detach a factor z. 

It follows from Eqs. (PT|) and that the restricted partition functions satisfy the boundary conditions: Zi{x = 
1, z, N) — and Zo{x — 0, z, N) — 0. For large N, one expects a separation of variables between x and N of the form, 

Zi,o(:E,z,iV)~A-^/i,o(x) (23) 

where A(z) is a function of z only (but independent of x) and the functions fi{x) and /o(x) satisfy the boundary 
conditions: /i(l) = and /o(0) = 0. Substituting this ansatz in Eqs. (PT|) and (1^^ and subsequently differentiating 
with respect to x, we get a pair of differential equations 

^ = -Xf,{x) - zXfoix) (24) 

^ = X.h{x) + Xfo{x) (25) 

Diagonalizing the [2 x 2] matrix, one obtains the solutions 

fi{x) = ae^^^^ + 6e-^^^^ (26) 

where a and b are arbitrary constants. The two boundary conditions /i(l) ~ and /o(0) = yield two relations 
between a and b 

^gAv^_^^g-AvT^ = (28) 
6^^^+-. ^ -e-^"^ = 0. (29) 



(1 - VT^) (1 + VT^) 

Eliminating a and b between Eqs. ([^5)1 and determines the eigenvalue A(z) exactly 



X{z) = / In 



1 + 



for < z < 1 (30) 



= --^i==tan"i [VI^T] for z > I. (31) 

Note that the function X{z) is analytic at z = 1. The form in Eq. (|3ip is just an analytical continuation of the form 
in Eq. jSO]) for z > 1. 

Substituting the large N form of Zi and Zq in Eq. (|20p . one obtains the large behavior of the partition function, 

Z(z, iV) = ^ P(M, N) z*^ - [A(z)]~^ (32) 

M 

where A(z) is given in Eqs. dSO]) and (l31|). Substituting the anticipated form of P{M,N) ^ exp[-A^<i>(A//, A^)] in Eq. 
(|32l) one gets 

Z(z, TV) = ^ P(M, iV) z^-f - f dyexp[^N {^{y)-y In z)]r^[X{z)]-^. (33) 

M 

Taking the N ^ oo limit in Eq. ((33|) gives 

min,, [$(y)-ylnz] = ln(A(z)). (34) 
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FIG. 1: The large deviation function ^{y) for < y < 1/2 in 1-d obtained from Eq. (I35p using Mathematica. 

Inverting via the Legendre transform finally gives the large deviation function 

= maxz [ln(A(z)) + y In z] , (35) 

where A(z) is given in Eqs. ([50)1 and ([51]) . Note that determining $(?/) from Eq. ([55)1 requires a knowledge of A(z) 
for all z > 0, thus we need both formulae of A(z) in Eqs. pO)) and ([3T|l . 

We have obtained $(?/) from Eq. (|35p using Mathematica and it is displayed in Fig.[T] Since the maximal value of 
M in 1-d is N/2, the allowed range of ?/ is < ?/ < 1/2. One can analytically obtain the form of <&(?/) in the three 
limiting cases y — > 0, y — > 1/3 and y 1/2. First let us consider the limit y 0. To find <&(?/) in this limit we 
need to use the z ^ form of A(z) in Eq. ([55)1 . As 2: ^ 0, it is easy to see from Eq. ([5(I|l that A(z) ^ ln(4/z)/2. 
Substituting this behavior in Eq. ([55|l and subsequently maximizing the r.h.s of Eq. ([55]) we find that <&(?/) diverges 
logarithmically, $(y) « In [ln2/(ej/)] as y — > 0. Next we consider the y — > 1/3 limit. Note that the mean number of 
minima (M) = 7V/3 in d = 1 which follows from the general result in Eq. ([7]). Thus y ^ 1/3 limit corresponds to 
behavior of M near its mean and one expects a quadratic form for ${y) near y = 1/3. Indeed this also follows from 
Eq. ([55]) . The limit y ^ 1/3 corresponds to using the z ^ 1 behavior of A(z) in Eq. ([55]) . Substituting z = 1 + e 
in Eq. ([5T|) and expanding in powers of e, one gets ln(A(z)) = — e/3 + 13e^/90 + 0{e^). Substituting this result on 
the r.h.s of Eq. ([35]) and maximizing one gets the expected quadratic behavior, <I>(y) « 45(y — 1/3)^/4 near y ^ 1/3. 
This is thus a special case in 1-d of the general behavior in Eq. ([H)l with a = 1/3 and b = 2/45. Finally, to derive 
the maximally packed limit y 1/2, we need to use the z — s- cx) behavior of A(z) in Eq. (|35p. In this limit, it follows 
from Eq. ([5T|l that ln(A(z)) ln(7r/2) — ln(z)/2 — 2z^^/^/t:. Using this form and maximizing the r.h.s. of Eq. ([55)1 
we get $(y) ~ ln(7r/2) -I- 2(1/2 — y) ln(7r(l/2 — y)/e) as y ^ 1/2. Thus, summarizing the three limiting behaviors 

(f>(y) w ln[ln2/(ey)] as y^O 
45 

« _(y- 1/3)2 as y-^1/3 

w ln(7r/2) + 2(1/2 -y) ln(7r(l/2-y)/e) as y^l/2. (36) 

Note that as one approaches the maximally packed limit y 1/2, ^(ymax = 1/2) = ln(7r/2). Thus, it follows from 
Eq. dni) that in 1-d, for large N 

-P(AW = A^/2,iV)~7-^, with 7 = V2, (37) 

the result declared in Eq. ([TB]) . 



IV. PROBABILITY OF THE MAXIMALLY PACKED CONFIGURATION: SOLVABLE CASES 

In this section, we focus only on the maximally packed configuration (where the number of minima on the lattice is 
maximal). From the general large deviation theory, we have already argued that the probability of such a configuration 
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P{Mynax, N) is expected to decay exponentially with the system size N as in Eq. ((T5)) . The goal is to compute the 
nontrivial constant 7. In the previous section, we have shown that in a 1-d chain, 7 = tt/2. In this section, we 
compute 7 exactly in a few other solvable cases, notably for a Cayley tree with fi branches and also for a two-leg 
ladder in 2-d. 



A. Exact Calculation of 7 on a Cayley Tree 

We consider a Cayley tree with /i branches and n generations. We label the generations by / = 1, 2, . . . , n starting 
from the leaf sites at the bottom (see Fig. [2]). The total number of sites on the tree is 

N = l + ^I + ^J.^ + ... + ^i"-^ = ^^—^. (38) 

Note that in the limit /k — > 1+, the tree reduces to a 1-d chain with N = n sites. 



1=4 



1=3 



1=2 



1=1 
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FIG. 2: The maximally packed configuration of local minima (denoted by black dots) on a Cayley tree with /i = 2 branches 
and n = 4 generations. The layers are labelled / = 1, 2 . . . starting with the bottom layer. 



At each site of this tree resides a random number Xi drawn independently from the uniform distribution over [0, 1]. 
We want to calculate the probability of the configuration with the maximum possible number of local minima. Note 
that if a site is a local minimum, its neighbours cannot be local minima. To find the maximally packed configuration, 
we note that the number of sites at the bottom layer {I = 1) increases exponentially with n as fi^~-^. Thus, to achieve 
maximal packing, it is necessary to fill up the bottom layer (the leaf sites) with local minima. Then the layer just 
above the bottom layer (/ = 2) is devoid of local minima. The next layer {I = 3) can again be packed with local 
minima. Thus, the maximally packed configuration is the one where alternate layers are fully packed with minima, 
starting with a fully packed layer at the bottom, as shown in Fig. [51 The total number of minima Mmax m this 
maximally packed configuration depends on whether the number of generations n is even or odd 



n+l 



1 



fl^ - 1 
/i2- 1 



even n 



odd 



(39) 
(40) 



For n even, the root is not a minimum whereas for n odd, the root is a minimum. In either case, for n large, the total 
number of minima in the maximally packed configuration is proportional to the total number of sites N , 



Mr, 



M+ 1 



-N. 



(41) 



Having identified the maximally packed configuration, we will next compute the probability of its occurrence and 
show that for large N 



where 



1 



7 = 



B 



1 



1 



/i + 1 ' /i + 1 



(42) 
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For simplicity, let us assume that n is even. One can perform an identical calculation when n is odd. For n even, 
all the even layers (labelled by I — 2m with m — 1,2, ... , n/2) are devoid of local minima whereas all the odd layers 
(labelled by / = 2m — 1 with m — 1,2,..., n/2) are fully packed with local minima. Our aim is to write a recursion 
relation. For this, it is convenient to define two probabilities P2m{x) and Q2m-i{x) defined respectively for even and 
odd layers. We define P2m{x) as the probability that a subtree of 2m generations (counted from the bottom of the 
tree) has a random variable x at its root and that x is not a local minimum. On the other hand, Q2'm-i{x) is the 
probability that a subtree of (2m — 1) generations (again counted from the bottom of the tree) has a random variable 
X at its root and that a; is a local minima. It is then easy to see that they satisfy the recursion relations 



P27n{x) 
Q2m+l{x) 



Q2m-iiy) dy 



P2m{y) dy 



(43) 
(44) 



where m = 1,2, ... , n/2. The recursions start with the initial condition Qi{x) = 1. Since n is even, the root of the 
full tree with n generations is not a local minima and hence the probability of the full tree is just P„(x) given that 
the value at the root is x. The probability of the maximally packed configuration is then obtained by integrating over 
X at the root 



P(Mmax, N) 



Pn{x)dx. 



(45) 



Thus, we need to solve Eqs. (|43|) and (|44|) and then substitute the solution for Pn{x) in Eq. (|45| to calculate 
P{M^,^,N). 

To solve the nonlinear recursion relations, it is convenient to define P2m{x) — [P2mix)]^^^ and q2m+i{x) = 
[Q2n.+i{x)]'^''. Then the recursions in Eqs. (|43p and (|44p become 



P2mix) 



(l2m-l{y) dy 



q2m+i{x) = / P2miy)dy 



(46) 
(47) 



starting with qi{x) = 1 for < a; < 1. It follows from Eqs. ([46]) and (|47|) that they satisfy the boundary conditions, 
P2m(0) = and (Z2m+i(l) — for all m > 1. However, P2m(l) and g2m+i(0) are nonzero. It is then useful to 
define the ratios, f2m{x) = P2m(a;)/p2m(l) and g2m+i{x) = q2m+i{x)/q2m+i{0) so that /2m(l) = 1, /2m(0) = and 
.92m4-i(0) = 1, g2m+i(l) = 0. In tcrms of the new functions, the recursions become 



f2m{x) 
g2in+l{x) 



P2m{l) JO 
'72m+l(0) 



52m- 1(2/) dy 



f2miy)dy 



(48) 
(49) 



starting with gi{x) — 1 for < a; < 1. As m increases, we expect that the ratio functions f2m{x) and g2m+i{x) will 
approach their respective m-independent fixed point forms f{x) and g{x). This means that as m - ~" ^2^-1(0) 



and 



p^,„(i) 

92,^+1(0) 



' P2r,i(l) 



A? and 



fix) = Ai / g'^iy)dy 
Jo 

9{x) = A2 / r{y)dy 



Ai 

(50) 
(51) 



with the boundary conditions /(O) = 0, /(I) ~ 1 and g{Q) — 1, g{l) — 0. These boundary conditions will determine 
the eigenvalues Ai and A2. 

Differentiating Eqs. (|50|) and (|5ip with respect to x gives 



dx 



dg 
dx 



-- Mg^ix) 

-\2F{X). 



(52) 
(53) 
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Multiplying ([52)1 by A2/^(x) and (|53p by Xig^{x), adding and then integrating, we find a conserved quantity 

X2r+\x) + X,g^^+\x) = C (54) 

where C is a constant independent of x. Putting a; = 0, 1 and using the respective boundary conditions gives 
C = Ai = A2. Thus, Ai = A2 = A and 



f+^x) + gf'+^x) = 1. 

The common eigenvalue A is yet to be determined. Eliminating g{x) between Eqs. ((52|) and ([55|) gives 

ax '- ■' 

subject to the boundary conditions /(O) = and /(I) = 1. Integrating and using /(O) = we get 

/•^(^) dz 

Jo ~ ~ 



(1 - 2M+1)m/(m+1) 

Using the other boundary condition /(I) — 1 determines A explicitly 



= Ax. 



A = 



B 



(55) 



(56) 



(57) 



(58) 



where B{m,n) = Jq x"^ ^(1 ^ 2;)" ^c?x is the standard Beta function. The eigenfunction f{x) in Eq. ([57]) can be 
expressed as the solution of the equation 



f{x)F 



fi+i /i+i fi+i 



B 



1 



1 



(59) 



where F[a, 6, c, z] is the hypergeometric function. The other eigenfunction g{x) follows from Eq. (j55p . or simply from 
the symmetry g{x) = f{l — x). To check that the solutions to the original recursion relations (|46p and P7|) indeed 
converge to these fixed solutions, we have numerically solved Eqs. and (|¥7)l for /i = 2. We find that the numerical 
solutions (the ratio functions f2m{x) and g2m+i{x)) converge to fixed point functions rather quickly after about 3 or 
4 iterations. In Fig. I52| we compare the numerical fixed point solution (the solution after 10 iterations) f{x) with the 
analytical solution in Eq. (|59p with fi ~ 2. The agreement is perfect. 



0.5 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 3: The numerically obtained fixed point function f{x) after 10 iterations, for fj, — 2, is compared to the analytical scaling 
function in Eq. (|59p obtained using Mathematica. The two curves are difficult to distinguish indicating perfect agreement. 



To determine 7, we use Eq. (|45p and the recursion relation (|47p which show that 



P(M^ax, N) 



Pn{x)dx = (J„+l(0). 



(60) 
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Then we take ratios to extract the bulk contribution of the Cayley tree, thereby removing artifacts coming from its 



surface [2J]: as m — > oo, we have 



P2m{l) 



A and 



g2m+l(0) 



+i(o) ^ ^' Ehminating P2m(l) gives a recursion for large m 
1 



r'?2m-l(0)- 



(61) 



Iterating it, we obtain, for large m, g2m+i(0) A ^)/(t^ i). Substituting this result in Eq. ([60|) and using 

N = (/i" — l)/(/i — 1) gives our final result for large N 



P{Mn,^^, TV) - 7 



-AT 



with 7 = A 



B 



(62) 



Note that the recursion relation (j6ip is valid only for large m. But in deriving the results above we have assumed 
that it holds even for small m. Even though the asymptotic result is not expected to change due to this 'initial 
condition' effect, one can avoid this "surface effect" by appropriately defining 7 as the ratio 



— = lim 

7 n — ^00 



Jo 




dx 




Jo Pn{x)dx 





(63) 



This definition follows from the following observation. The numerator on the r.h.s of (f63|) denotes the probability of 
a maximally packed configuration on a tree with (n + 1) generations. The denominator is the joint probability that 
the disconnected subtrees of n generations are all maximally packed. Thus, the ratio on the r.h.s. is just the factor 
by which the probability of a maximally packed configuration per site decreases when one fuses the /i number of n- 
generation trees and one additional root to construct a newly maximally packed configuration on a (n+ l)-generation 
tree. But, asymptotically for large n, this is precisely I/7 by the orginal definition in Eq. p^ . Physically, ln7 is the 
additional entropy change when packing an extra minimum in the tree. Using the recursion relations (|43p . (j44p in 
Eq. one gets 



— — lim 

7 n — >oo ^ 



/■ 




dx = 




.'7«+i(0). 


Jo 



9n+ii.x)dx. 



Substituting g^{x) — 4 df /dx, integrating and using /(I) 



7 



1 



B 



1, we get 
1 



1 



+ 1 > + 1 



(64) 



(65) 



Note that for = 1, we recover the 1-d result, 7 — 7r/2. For ^ = 2, we get 7 = 5(1/3, l/3)/3 = 1.76664 .... As 
/i — > cx), 7 converges slowly to 7 = 2. Since the /i — > cx) result should coincide with that on a hypercubic lattice in 
the infinite dimension limit, we expect that for hypercubic lattices in d dimensions, 7 is bounded, tt/2 < 7 < 2. As 
d increases from 1 to 00, 7 should increase monotonically from its d = 1 value 7r/2 = 1.56079 ... to 2. Our numerical 
solutions on d-dimensional hypercubic lattices with d = 2, 3, 4, 5 are consistent with these bounds. 



B. Maximally Packed Configuration on a Bipartite lattice: An Equivalent Plaquette Model 

On a bipartite lattice, the maximally packed configuration is the one where one places a local minimum at every 
alternate site. For example, on a square lattice, a maximally packed configuration has a checkerboard pattern as 
shown in Fig. |4l Let us first focus on the neighbourhood of a single local minimum. The sites at the corners of 
the plaquette containing this minimum are clearly not local minima. Let 0:1, X2, x^ and x^ denote the values of the 
random variables at the four corners of the plaquette. Then, given these four values, the probability that the site at 
the center of the plaquette is a local minimum is clearly 

Pi{xi,X2,X3,X4) = minixi,X2,X3,X4:) (66) 

where we have used the fact that the random variable at the center of the plaquette is drawn from a uniform 
distribution. So, the probability of the full checkerboard configuration where every alternate site is a local minimum. 
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given the values of the random variables {x^j's at the sites of the other sublattice, is obtained by multiplying all the 
plaquettes 

Prob (M^ax, N I {x,}) - II min {xi{j),X2{j), X3{j),X4{j)) (67) 

plaqucttc J 

where the plaquettes are labelled by j and xi{j), X2{j), X3{j) and X4{j) are the four random variables at the corners 
of the j-th plaquette. The total number of plaquettes will also be denoted by M since the number of plaquettes is the 
same as the number of local minima. To simplify, it is convenient to rotate the lattice by —45° (as shown in Fig.[4|). 



-45° rotation 




-<) — • — Q- 



Q — » — ,a- 



-X) 



-a — V 



a — • — — f — 6 



-O o o 



6 6 9 ^ 



-6 Q o <? 



6 G 6 o 



FIG. 4: A maximally packed configuration on a square lattice. On this checkerboard pattern, the black circles denote the 
location of the local minima. The corners of every plaquette around a minima cannot contain a minimum and are shown by 
empty circles. On the right, the same pattern, rotated by an angle —45°. 



We then have a plaquette model where at each corner of a plaquette lives a random variable Xi drawn from a 
uniform distribution over [0, 1] and we are interested in calculating the product in Eq. (|67p . Finally, the probability 
of the maximally packed configuration on the original lattice is obtained by averaging over the x^'s with uniform 
distribution in [0, 1] 

P{M^.^,N) = { II mm{xiij),X2{j),X3{j),Xi{j))) (68) 

plaquette j 

where the angled brackets simply indicates integration over all the x variables from to 1. Note that this plaquette 
model is very general and can be extended to any bipartite lattice. Also, the right hand side of Eq. (|68p is actually the 
probability of just one of the two checkerboard configurations, so the true -P(Af„iax: N) is actually twice this amount, 
but to avoid unnecessary complications we shall keep to this notation. 



1. Plaquette Model in One Dimension 



As a simple example, for a 1-d chain, one can reproduce the result 7 = 7r/2 quite easily using the plaquette 
representation. The number of plaquettes is clearly M — N/2. In this case, Eq. ([55)1 gives 

M=N/2 

mm{xi,Xi+i)) (69) 

i=l 

where i runs over every alternate sites of the original 1-d chain. The quantity on the r.h.s. of (|69p can be evaluated 
using a simple transfer matrix approach. Defining a transfer matrix via, {xi\T\xi+i) = min{xi, Xi^i), we have from 
Eq. ([69|) . assuming a closed periodic chain 



P(M„,ax, N) = Tr [f (70) 
where Tr is the trace. The eigenvalue equation of the transfer matrix is 

r-l 

mm{x,y)tp{y)dy ^ Xipix) (71) 
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where ijj{x) is an eigenfunction with eigenvalue A. Dividing the range of integration into [0,2:] and [x, 1] and differen- 
tiating twice, one gets 



^S + ^(-) = (72) 



with the boundary conditions ip{x = 0) = and = 1) = where ip'{x) = dip/dx. The solution is simply, 

■0(2;) = Asm{x/^/X) where 

^- .^(2m\l)^ ^ --O^M... (73) 
The largest eigenvalue is A = A/tt'^ corresponding to m = 0. Thus for large M — N/2, one gets from Eq. ([70)1 



P(M,„ax,iV)^ 

reproducing the results obtained in the previous section. 



N/2 

- k/2]-^ (74) 



2. Plaquette Model on a Two-leg Ladder 

Another nontrivial solvable case of the plaquette model is on a two-leg ladder shown in Fig. [5l It can be thought 
of as the first layer of the full 2-d model shown on the right in Fig. [H At the center of each plaquette lives a local 
minimum. Let M be the number of plaquettes. Then the total number of sites in the lattice (counting the minima at 
the centers) is = 3M. 




TWO-LEG LADDER 



FIG. 5: The plaquette model on a two-leg ladder. At the center of each plaquette lives a local minimum shown by the black 
circles. 

The probability of the maximally packed configuration in Eq. (j67p can again be computed by defining the transfer 
matrix, {xi, X2\T\x3, X4) = min(a;i, 2:2, 2:3, 2:4) where (2:1,2:2) refers to the two corners on the left of a plaquette and 
(2:3,2:4) refers to the two corners on the right of a plaquette. Then, 

P(M^ax, A^) = Tr [f (75) 

The corresponding eigenvalue equation is now two dimensional 

1 .1 

o?2;3 / 1^2:4 min(2;i, 2:2, 2:3, 2:4) -0(2:3. 2:4) = A ■0(2:1, 2:2) (76) 
^0 

which is considerably harder to solve compared to the 1-d case. While one can solve the 2-d eigenvalue equation 
directly, it is somewhat easier to first reduce it to an equivalent 1-d eigenvalue problem by using the following trick. 
Let us first define the random variable Zi = min(2:i (i), 2;2(«), 2:3(2), 2:4(1)) at the center of each plaquette i. Clearly 

P(Afmax,A) = (ziZ2...Zm) (77) 
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FIG. 6: On the left we have a single plaquette with the minimum zi — min(a;i, 3:2, 3:3, 2:4) at the center. On the right we 
have two adjacent plaquettes with minima zi — min(a;i, 3:2, 3:3, 2:4) and Z2 = min(3;3, 3:4, 3:5, 3;6) at their respective centers. The 
variables 21 and Z2 are corrrelated as they share the common bond with the two elements 3:3 and 3:4. 



where M = N/3 is the number of plaquettes. The random variables ZiS are obviously correlated as two adjacent 
plaquettes will share two common random variables as shown in Fig. [6l 

To evaluate the average on the r.h.s of Eq. (|77p we need to know the joint distribution of the ^^'s which are 
correlated (see Fig. This joint distribution can be explicitly computed. To see this, let us first compute the 
cumulative distribution of one single z variable, say zi = min(xi(l), 0:2(1), 2:3(1), a;4(l)). For simplicity, we denote 
xi(l) = xi, 2:2(1) — X2 etc. (see Fig. [6]). Evidently 

Prob(zi > y) ~ Prob(min(2:i, X2, 2:3, 2:4) > y) = (1 — y)* (78) 

where we have used the fact that xi, X2, 2:3 and 2:4 are all independent random variables drawn from the uniform 
distribution over [0, 1]. Next, let us consider the joint distribution of two consecutive z^s, say zi and Z2- Let us denote 
the random variables at the corners by 2:1, 2:2, X3, 2:4, 2:5 and xg as shown in Fig. [51 Then it is easy to see that 

Prob(zi > yi, Z2 > V2) = Prob(min(xi, 2:2, 2:3, 2:4) > j/i, min(x3, 2:4, 2:5, 2:6) > ^2) 

= (1 - [1 - max(yi, y2)]' (1 - 2/2)'- (79) 

In a similar way, one can construct the joint distribution of three adjacent plaquette minima 

Prob(zi > yi, Z2 > y2, 23 > ^3) = (1 - Vif [1 - max(yi, y2)]^ [1 - max(y2, y3)]^ (1 - y3)^- (80) 

One can repeat the process above for higher number of adjacent plaquettes and one immediately sees the pattern for 
the full ladder. Let us denote 

F{yi,y2,y3, ■ ■ -.Vm) = Prob(2;i > yi, Z2 > y2, ^3 > ya, ■ ■ ■ , zm > yn)- (81) 
Then, for a ladder with open ends, we get 

■ M 



^(yi,y2,y3, ■ • • ,2/m) (i - yif 



W[l max(yi_i,yj)]^ 



.1=2 



{^-VM?. (82) 



For a ladder with closed (periodic) ends, this is even simpler 



M 



F{yi,y2,y3, ■ ■ -.Vm) = H ~ max(yj_i,y,)]^ . (83) 



i=l 



where one identifies yo ~ vm- 

Once we have the joint distribution, it is easy to rewrite the average on the r.h.s. of Eq. (j77p in terms of the joint 
distribution 



P{M^^^,N) ^ {ziZ2. ■ -zm) ^ j .../ dyidy2 . ■ -dyM F{yi,y2,y3, 

Jo Jo 



,yM) (84) 
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The latter identity can be easily derived using integration by parts. Let us, for simplicity, consider a plaquette with 
closed ends. Then, substituting the joint distribution from Eq. ([55)1 into Eq. we now have a one-dimensional 

multiple integral to peform 

P{Mn,,^, N)= ... dyidy2 . . . dyM TT " max(y,_i, y,)]' . (85) 
Jo Jo 

Using the identity 1 — max(?/i, 1/2) = min(l — yi, 1 — 2/2) and making a change of variable Xi ^ I ~ yi, the integral in 
Eq. ([85]) simplifies further 

P{M^a,^, N) ^ ... j dxidx2 . . . dxM TT [min^(xi_i, x^)] . (86) 
^0 Jo 

This one dimensional integral can now be performed by a transfer matrix technique which operates in one dimension. 



Defining {x.i^i\T\x^) = min^(a;^_i , a;^) we get, P{M^^^,N) = Tr 



T 



M 



and the eigenvalue equation is given by 



rnin^ {x, y) 'ilj{y) dy = XiIj{x). (87) 

Thus we have managed to reduce a two-dimensional eigenvalue problem in Eq. (|76p into an equivalent but much 
simpler 1-d eigenvalue problem in Eq. (j87p which can then be solved exactly. 

To proceed, we first divide the range of integration in Eq. l(87)) into [0, x] and [x, 1]. Next we differentiate once to 
get 



Note the boundary conditions that emerge from Eas. ([57|) and 'tlj{x = 0) = and tp^x = 0) = 0. But in addition 
one also has to satisfy ip'{x — 1) — 0. Dividing Eq. ((88|) by x and differentiating once more we get an ordinary second 
order differential equation 

^"{x) - -^P'ix) + 4- VX^) = (89) 

X A 

with the boundary conditions: ^p{x = 0) = 0, ^p'{x = 0) = and ^p'{x = 1) = 0. The general solution to this equation, 
after a few changes of variables, can fortunately be obtained explicitly as a linear combination of two independent 
Bessel functions. The boundary condition ip{0) = rules one of them out. Then, the most general solution satisfying 
-0(0) = can be written as 



8x3 
9A 



^ljix)=AxJ2/3\\l—] (90) 



where Ju{x) is the ordinary Bessel function with index v ^25] and A is an arbitrary amplitude. Note that the other 
boundary condition at a; = 0, namely V''(0) = is automatically satisfied by the solution in Eq. ((90)) . This can be 
seen by using the small x expansion of J,y{x) ~ x'^ which indicates ^/'(cc) ~ a;^ as a; ^ 0. Hence V''(0) = 0. 

To determine the eigenvalue A, we have to use the other nontrivial boundary condition at x = 1, namely ip^x — 
1) = 0. This condition, substituted in Eq. ([50)1 . gives us an implicit equation for A that looks a bit complicated 

However, a nice simplication occurs when one uses the identity [1^, xJl(x) + vJi,{x) = a;J^_i(a;). Then Eq. (pij) 
simply gives 
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The Bessel function oscillates on the positive axis of its argument, so each of its zeroes (roots) would give an eigenvalue 
A. However, for large M , we are interested only in the largest eigenvalue which is then given by 



A 



9a2 



(93) 



where a is smallest positive root of J-i/3(x) = 0. The root a is known 26], a — 1.86635 . 

Using P(M,„ax, N) = Tr 
ladder 



T 



M 



A for large M and using M — N/i, we finally get the exact result for the two-leg 



P(M„ax, N) ^ J 



-N. 



where 7 



9a^ 



1/3 



1.57657. 



(94) 



When compared with the 1-d result, 7 = tt/2 = 1.57079 . . ., we see that 7 changes by a very small amount as one 
goes from a chain to a ladder. It would be interesting to see if one can extend these calculations to ladders with more 
than two legs and eventually to the full 2-dimensional lattice. 



V. NUMERICAL RESULTS FOR MAXIMALLY PACKED d-DIMENSIONAL LATTICES 

Now we consider the dependence of 7 on dimension. For simplicity, we have focused on d-dimensional hypercubic 
lattices as these lattices are simple to parametrize and are bipartite, allowing us to use the d-dimensional analog of 
the framework described in Sect. IIVBI Our approach will be computational: we numerically estimate -P(Mmax, N) 
from an integral representation for lattices of increasing size and then try to extract the limit of large sizes. To avoid 
boundary effects which slow down the convergence to this limit, we have used L'^ lattices having periodic boundary 
conditions in all directions. 

The integral in Eq. (|68p is based on decomposing the lattice into even and odd sites and imposing the local minima 
to be say on the even sites. Its generalization to d dimensions involves all the neighbors of a given even site j, the 
corresponding "star" set of 2d odd sites playing the role the plaquette had for the square lattice: 

P(M„,ax, ^) = ( n "™ • ■ • X2dij)) ) . (95) 

star j 

Here the average is over all values of the random variables belonging to stars; these are uniform i.i.d. in [0,1]. 
Furthermore, as before, we neglect the factor 2 in /'(A/max, N) coming from the fact that the minima could have been 
taken to be on the odd sites. The difficulty in computing these integrals is their high dimensionality. In Monte Carlo 
as used in most statistical physics applications, it is straightforward to use importance sampling methods (e.g., the 
Metropolis algorithm) to get expectation values of observables; unfortunately here, the quantity to compute is the 
analog of the free energy and it is not directly accessible via such methods. We have thus used a different approach 
that is of the "variance reduction" type. It can be motivated as follows. In Eq. ((95)) we are to get the mean value of 
the integrand, sampling the random variables Xi uniformly. It is quite easy to see that the signal to noise ratio goes 
to zero exponentially with the number of lattice sites; to counteract this, we sample the Xi with a different density 
and then correct for this biased sampling. For this to be practical, we keep the Xi as i.i.d. variables, but optimize 
their individual distribution so as to maximize the signal to noise ratio. Let p{x) be the probability density used for 
sampling the Xi. For any such distribution. 



P(M_, A^) = ( n "^^^(-^(^•)'^^(^-)--^^(^-)) ) (96) 

starj [piMj))pix2{j)) ■ ■ ■ P{x2dij))] 



The denominator corrects for the modified measure of the random variables, taking into account the fact that each 
odd site appears in 2d stars. If p is well chosen, the numerator and denominator of the integrand will fluctuate 
together so that their ratio has a reduced variance. For our purposes, we parametrized p as follows: 

p{x) = Ax'^ X < X* 

p{x) = B x>x* (97) 

where (j) and x* are arbitrary parameters while A and B are set so that p is continuous and is a normalized probability 
density. For each dimension, we adjust (j) and x* to minimize the variance the integrand; the signal to noise ratio 
still decreases exponentially with the number of lattice sites, but with a smaller exponent. The motivation for this 
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d=l d=2 d=3 d=4 d=5 

-yhypercube 1.57082(6) 1.6577(6) 1.7152(10) 1.761(5) 1.806(10) 
-ycayiey 1.570796 1.854075 1.927621 1.956922 1.971464 

TABLE L The numerical estimates of 7 for d-dimensional hypercubic lattices {'yhypercube)', error estimates come from statistical 
fluctuations and uncertainties in the large lattice size extrapolation. Bottom line: exact values for the Cayley tree with the 
same coordination number, n + 1 = 2d. 

functional form is simple: the odd sites, all of which are maxima, have an a posteriori distribution that is strongly 
suppressed at low values of the random variable. 

One last obstacle comes from the fact that the integrand typically takes on small values when one takes large lattices; 
this is expected of course since the integral itself is becoming exponentially small. To keep track of values that are 
smaller than what can be represented by the machine coding, (we used 96 bit representations of real numbers), we 
shifted multiplicatively each star term in the integrand and corrected for this shift when computing ln(P(Mmax, .^)). 
The set of these procedures then gave us values of ln(P(Mmax, N))/N with measurable statistical errors for a range 
oi N = L'^ from which we extrapolated to the large N limit. 

In practice, we find that this strategy works very well in low dimensions. For instance in dimension d = 1, we still 
have a very good precision for i = 30 sites, and the large L limit can be very reliably extracted, giving 7 to better 
than 5 significant figures. (Of course, since the exact value is known at d = 1, this really only serves as a check of 
our procedures.) In dimension d — 2, the method gives better than 4 significant figures for L w 16. Our estimates 
of 7(iV) = P{Mma.^, N)~^^^ as a function of 1/L — 1/Vn are shown in Fig. ([7]). The convergence to the L = 00 
limit seems to follow a 1/L law, a property we also find for the higher dimensions investigated. Here we see that 

1.68 



1.675 
1.67 



1.665 

0.05 0.1 0.15 0.2 

1/L 

FIG. 7: The convergence to the large L limit of the estimates of 7(A'') for the square lattice, N = . 

7(iV) converges rapidly to a value close to 1.658. Unfortunately the signal to noise ratio decreases with L and with 
dimension; furthermore the accessible range of L decreases fast as d increases; because of this, we are able to extract 
7 reliably only for dimensions up to d = 5. Our results are summarized in table U where the error estimates come 
from both statistical noise and uncertainties in the large L extrapolation. 

VI. DISCUSSION AND CONCLUSIONS 

We have been concerned with the statistical properties of M, the number of local minima in a random energy 
landscape. The low order moments of M can easily be obtained; higher order moments could be computed by 
automated counting of graphs. The atypical values of M follow a large deviation principle as given in Eq. I|13p : we 
have been able to compute this function $(j/) in one dimension analytically; it diverges as y — > and goes to a finite 
limit when M reaches its maximum value where one half of the lattice sites are local minima. This then led us to 



17 



consider the limit of maximum packing (M takes on its largest possible value) for more general lattices. We focused 
on bipartite lattices where up to half of the sites can be local minima. We derived analytically the probability of 
this maximum packing P(Mi„ax, N) for the Cayley tree and for a two-leg ladder. We then tackled d-dimensional 
hypercubic lattices by computational techniques. For all these lattices, it is easy to see that P(Mmax, N) > 
simply by forcing the Xi on the even (odd) sites to be less (greater) than 1/2; this immediately leads to 7(d) < 2 for 
all dimension. Furthermore when d becomes large, it will be very rare to have maximum packing if the Xi do not very 
nearly satisfy this even-odd pattern so one expects ^{d) to tend towards 2 in the large d limit. 

Given this large d limit, It is natural to ask how 7 = 2 is approached. In Fig. ([S]) we show that the Cayley tree 
case follows a very clear power correction law which can be derived analytically as being 1/d^. The case of the 
d-dimensional lattice is less clear but is compatible with a 1/c? law. 

1 



0.1 



0.01 

0.1 1 

lid 

FIG. 8: The convergence to the large d limit of the estimates of 7(d) vs plotted is the difference 2 — 7(d) which vanishes 
as a power at large d. 

The Cayley tree thus provides an approximation to the hypercubic case but not a very accurate one; this is probably 
because the nature of the correlations from site to site of the Xi (given that one has a maximally packed configuration) 
are quite different when considering the tree rather than the hypercube. 

Our work can be extended in several ways, (i) We computed the exact large deviation function in one dimension, 
but this function can also be determined for the Cayley tree [2^. (ii) One can also introduce a chemical potential z 
for each minimum and consider the thermodynamics of minima as a function of z; recent numerical results by Derrida 
in 2-d 29] indicate the presence of a phase transition with z. (iii) Our variables Xi were i.i.d. random variables; does 
a large deviation principle still hold if these variables have short range correlations? One expects so. (iv) How do all 
these statistical properties generalize if one asks for minima within a given energy range? Finally, it would also be 
of interest to understand the similarities and differences between the statistics of local minima in continuous random 
energy landscapes [1, [s^l and in lattice models as presented here; the energies could be random or correlated as in 
the Sinai problem. 
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